Stochastic Gravitational Wave Background from Light Cosmic Strings 
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Spectra of the stochastic gravitational wave backgrounds from cosmic strings are calculated and 
compared with present and future experimental limits. Motivated by theoretical expectations of 
light cosmic strings in superstring cosmology, improvements in experimental sensitivity, and recent 
demonstrations of large, stable loop formation from a primordial network, this study explores a 
new range of string parameters with masses lighter than previously investigated. A standard "one- 
scale" model for string loop formation is assumed. Background spectra are calculated numerically for 
dimensionless string tensions G^/(? between 10~^ and 10~^*, and initial loop sizes as a fraction of the 
Hubble radius a from 0.1 to 10~®. The spectra show a low frequency power-law tail, a broad spectral 
peak due to loops decaying at the present epoch (including frequencies higher than their fundamental 
mode, and radiation associated with cusps), and a flat (constant energy density) spectrum at high 
frequencies due to radiation from loops that decayed during the radiation-dominated era. The string 
spectrum is distinctive and unlike any other known source. The peak of the spectrum for light strings 
appears at high frequencies, significantly afi^ecting predicted signals. The spectra of the cosmic string 
backgrounds are compared with current millisecond pulsar limits and Laser Interferometer Space 
Antenna (LISA) sensitivity curves. For models with large stable loops [a = 0.1), current pulsar- 
timing limits exclude Gptjc? > 10~^, a much tighter limit on string tension than achievable with 
other techniques, and within the range of current models based on brane inflation. LISA may detect 



a background from strings as light as G^/(? 
at roughly lO" GeV. 



10 , corresponding to field-theory strings formed 



PACS numbers: 11.27.-|-d, 98.80.Cq, 98.70.Vc, 04.30. Db 



I. INTRODUCTION 



It is well established in gauge field theories that the 
vacuum expectation value of a field {(j)) can take on non- 
zero values when the ground state, at the minimum of 
the zero-temperature potential or "true vacuum" , breaks 
the symmetry of the underlying Lagrangian P, In 
the hot early universe, thermal effects lead to an ef- 
fective potential with a temperature dependence, such 
that the vacuum starts in a symmetric "false" vacuum 
and only later cools to its ground state below a cer- 
tain critical temperature [11, d, 0, H, @| ■ In general how- 
ever the expansion of the universe occurs too quickly for 
the system to find its true ground state at all points in 
space Depending upon the topology of the mani- 
fold of degenerate vacua, topologically stable defects form 
in this process, such as domain walls, cosmic strings, 
and monopoles [H, H, S Hi- In particular, any field 
theory with a broken U(l) symmetry will have classi- 
cal solutions extended in one dimension, and in cosmol- 
ogy these structures generically form a cosmic network 
of macroscopic, quasi-stable strings that steadily unrav- 
els but survives to the present day, losing energy primar- 
ily by gravitational radiation E [III E Hill- A 
dual superstring description for this physics is given in 
terms of one-dimensional branes, such as D-strings and 
F-strings [H, [H, [B [H, [H H . 
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In this paper we calculate the spectrum of background 
gravitational radiation from cosmologies with strings. 
We undertake calculations at much lower string masses 
than previously, with several motivations: 

1. Recent advances in millisecond pulsar timing have 
reached new levels of precision and are providing 
better limits on low frequency backgrounds; the cal- 
culations presented here provide a precise connec- 
tion between the background limits a,nd fundamen- 
tal theories of strings and inflation |2l|, [^, 0, 0] . 

2. The Laser Interferometer Space Antenna (LISA) 
will provide much more sensitive limits over a broad 
band around millihertz frequencies. Calculations 
have not previously been made for theories in this 
band of sensitivity and frequency, and are needed 
since the background sp ectrum depends signifi- 
cantly on string mass pH . l26l . Wi\ . 

3. Recent studies of string network behavior strongly 
suggest (though they have not yet proven defini- 
tively) a high rate of formation of stable string 
loops comparable to the size of the cosmic horizon. 
This results in a higher net production of gravita- 
tional waves since the loops of a given size, forming 
earlier, have a higher space density. The more in- 
tense background means experiments are sensitive 
to lower string masses [H, [H, [s^l • 

4. Extending the observational probes to the light 
string regime is an important constraint on field 
theories and superstring cosmology far below the 
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Planck scale. The current calculation provides 
a quantitative bridge between the parameters of 
the fundamental theory (especially, the string ten- 
sion), and the properties of the observable back- 
ground gili. 

The results of our study confirm the estimates made in an 
earlier exploratory analysis, that also gives simple scaling 
laws, additional context and background [33|. Here we 
add a realistic concordance cosmology and detailed nu- 
merical integration of various effects. We estimate that 
the approximations used to derive the current results are 
reliable to better than about fifty percent, and are there- 
fore useful for comparison of fundamental theory with 
real data over a wide range of string parameters. 

Cosmic strings' astrophysical properties are strongly 
dependent upon two parameters: the dimensionless 
string tension G/i (in Planck units where c — 1 and 
G = mpi~'^) or GjjL/c^ in SI units, and the interchange 
probability p. Our main conclusion is that the current 
pulsar data 0, [H, [s^l already place far tighter con- 
straints on string tension than other arguments, such as 
microwave background anisotropy or gravitational lens- 
ing. Recent observations from WMAP and SDSS fill 
have put the value of G/x < 3.5 x 10~^ due to lack of cos- 
mic background anisotropy or structure formation con- 
sistent with heavier cosmic strings. From [sF] it is found 
that up to a maximum of 7% (at the 68% confidence 
level) of the microwave anisotropy can be cosmic strings. 
In other words, for strings light enough to be consistent 
with current pulsar limits, there is no observable effect 
other than their gravitational waves. The limit is already 
a powerful constraint on superstring and field theory cos- 
mologies. In the future, LISA will improve this limit by 
many orders of magnitude. 

In general, lower limits on the string tension constrain 
theories farther below the Planck scale. In field theories 
predicting string formation the string tension is related 
to the energy scale of the theory through the relation 
G^i cx A'^/nipi [H, 0, Hi- Our limits put an upper limit 
on As in Planck masses given by As < 10~'*'^, already 
in the regime associated with Grand Unification; future 
sensitivity from LISA will reach A^ 10^^, a range of- 
ten associated with a Peccei-Quinn scale, inflationary re- 
heating, or supersymmetric B-L breaking scales [36]. In 
the dual superstring view, some current brane cosmolo- 
gies predict that the string tension will lie in the range 
10"*^ < G^i < 10-" [211,121113; our limits are already 
in the predicted range and LISA's sensitivity will reach 
beyond their lower bound. 

Cosmic strings radiate the bulk of their gravitational 
radiation from loops (T], [ll|, [l^ |3^ . For a large network 
of loops we can ignore directionality for the bulk of their 
radiation which then forms a stochastic background [s^ . 
Each loop contributes power to the background: 



where 7 is generally given to be on the order of 50 to 100 

0, [13, EE El, m. 

Gravitational radiation from strings is characterized 
for our purposes by the dimensionless luminosity of grav- 
itational waves 7, which depends on the excitation spec- 
trum of typical radiating modes on string loops. Loop 
formation is characterized by a, the typical size of newly 
formed loops as a fraction of the Hubble scale H{t)^^. 
Estimates based on numerical simulations have in the 
past suggested very small values, leading to the hypoth- 
esis that a ~ G/i or smaller [isl . H . [45l |: those ideas 
suggest that all but a fraction ^ G^ or smaller of loops 
shatter into tiny pieces. Recent simulations designed to 
avoid numerical artifacts suggest a radically different pic- 
ture, that in fact a significant fraction of loops land on 
stable non-self-intersecting orbits at a fraction ^ 0.1 of 
the horizon scale [1^ [IE [13j [1E| • Our study is oriented 
towards this view, which leads to a larger density of loops 
and a more intense background for a given string mass. 

For units of the gravitational energy in a stochastic 
background a convenient measure we adopt is the con- 
ventional dimensionless quantity given by 
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where pgw is the energy density in the gravitational waves 
and Pc is the critical density. Q.gy^{f) gives the energy 
density of gravitational waves per log of the frequency in 
units of the critical density, is proportional to the 

mean square of the strain on spacetime from a stochastic 
background of gravitational waves. 

There has been much interest in cusps and kinks as 
sources of gravitational radiation from cosmic strings 
[47I , m, i^, ^] . In this paper an estimate of the radi- 
ation power from cusps and kinks is included as higher 
frequency behavior in the loop gravitational wave spec- 
trum, but the beaming and bursts of gravitational waves 
is not discussed explicitly. As estimated in [3^ . in the 
regime of very light strings, observable bursts are ex- 
pected to be rare, and harder to detect than the stochas- 
tic background. 

Finally we would like to state some additional uncer- 
tainties associated with the model of strings described 
below: 

• There is the possibility of multiple scales in loop 
formation: some string theories suggest multiple 
stable intersection points and loops connected by 
these intersections. This may strongly affect the 
accuracy of the one scale model, and indeed some 
of these theories also have significant other modes 
of decay that invalidate the entire framework used 
here (B!]. 

• The reconnection probability p can be less than 
one. This may affect loop formation and the 
gravitational wave background. In general it is 
thought that this effect always increases the the 
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predicted background although it is not agreed by 
how much p^ . 

• There may be transient behavior of cosmic strings 
due to the formation process of cosmic quench- 
ing. The model below assumes that transient ef- 
fects have long since settled out by the time our 
model begins; this seems likely at LISA frequencies 
which correspond to initial loop sizes much larger 
than the horizon at the end of inflation [18]. 

Further numerical string network studies will be needed 
to resolve these issues. 

II. MODEL OF STRING LOOP POPULATIONS 

A. Loop Formation 

The behavior of strings on cosmological scales has been 
thoroughly discussed in the literature 
0; Ea, nil- By the Kibble mechanism strings are 
created as a random walk with small coherence length, 
and they evolve following the Nambu-Goto action. The 
strings form in large threadlike networks [46j, and on 
scales longer than the Hubble length cH{t)~^ the net- 
work is frozen and stretches with expansion. The strings 
move with speed c, and they interact to form loops with 
a probability p by different mechanisms: one involves 
two strings, the other a single string forming a loop. 
These loops break off from the infinite string network 
and become the dominant source of gravitational waves. 
They oscillate with a fundamental frequency given by 
/ = 2c/ L, where L is the length of the loop. In general 
the loops form a discrete set of frequencies with higher 
modes given by /„ = 2nc/L. 

The exchange probability p measures the likelihood 
that two crossing string segments interact and form new 
connections to one another. If p=0 then the string seg- 
ments simply pass through each other; but if p=l then 
they can exchange "partners" , possibly forming a loop. 
The value of p depends upon the model [1^ Hj, [s^ ; it 
is close to unity in models most commonly discussed but 
in principle is an independent parameter (and in a broad 
class of models, the only one) that can modulate the am- 
plitude of the spectrum for a given string tension. In this 
paper it is taken to be p=l and is parameterized as part 
of the number of loops formed at a given time 0, [lo|,]33|. 
The number of loops is normalized to the previous results 
of R.R. Caldwell and B. Allen [s^, which is described in 
detail later in this section. 

B. Loop Radiation 

For the dynamics of a string the Nambu-Goto action 
is used: 

S^-H j ^f^dTda, (3) 



where r and a are parameterized coordinates on the 
world sheet of the loop, and g^'^^ is the determinant 
of the induced metric on the worldsheet. The energy- 
momentum tensor of the string can be given by [l| , 

T^''(t,x)=/i j da (^i^x" -^^^S^^H^-^it,a)), 

(4) 

where r is taken along the time direction. The radiated 
power from a loop at a particular frequency, /„, into a 
solid angle 17 is, 

^ = 4nGf^ (T;Mn,k)T'^-ifn,k) - \ |T;(/„,k)|^) . 

(5) 

Here T^^(/„,k) is the Fourier transform of Ttji/(i,x) 
Dimensionally the power is given by G^j? 13] with nu- 
merical coefficients determined by the frequency mode 
and the luminosity of the loop. 

In general for the power output of a loop at a given 
mode we can write En — PnG\j?c P, [s^, where di- 
rectionality is ignored by effectively inserting the loop 
into a large network of loops and the P„ are dimen- 
sionless power coefficients. Note that c has been rein- 
serted. The power coefficients are generally given the 
form Pn oc n"'*/'^ [1] which will create an /"^/'^ high fre- 
quency tail to the power spectrum of a loop. The total 
power from a loop is a sum of all the modes: 

oo 

E = ^P„Gm'c, (6a) 

ri=l 

= 7Gm^c. (6b) 

Thus 7 measures the output of energy and is labeled the 
"dimensionless gravitational wave luminosity" of a loop. 
For most calculations we use the value 7 = 50 [42] . 

For the majority of our calculations we assume only the 
fundamental mode contributes, effectively a delta func- 
tion for P{f), and / — 2c/L is the only emission fre- 
quency of each loop. This assumption is relatively good. 
To check this, some calculations include higher frequency 
contributions and show that the effect on the background 
is small overall, with almost no difference at high frequen- 
cies. This is discussed in more detail in Results. The 
high frequency additions also represent the contribution 
from kinks and cusps, which add to the high frequency 
region with an f~^l^ tail from each loop. This high fre- 
quency dependence is the same as that for the high fre- 
quency modes, so their contributions should be similar, 
particularly for light strings. Light strings contribute at 
higher frequencies because the loops decay more slowly 
and more small loops are able to contribute at their fun- 
damental frequency longer throughout cosmic history. 

The populations of loops contributing to the back- 
ground fall naturally into two categories [s^: the high 
redshift "H" population of long-decayed loops, and the 
present-day "P" population, where the redshift is of the 
order unity or less. The spectrum from the H population. 
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which have long decayed and whose radiation is highly 
redshifted, is nearly independent of frequency for an early 
universe with a p = p/3 equation of state. The P loops 
have a spectrum of sizes up to the Hubble length; they 
create a peak in the spectrum corresponding approxi- 
mately to the fundamental mode of the loops whose de- 
cay time is the Hubble time today, and a high frequency 
tail mainly from higher modes of those loops. For all of 
the observable loops (that is, less than tens of light years 
across), formation occurred during radiation era (that 
is, when the universe was less than 10^ years old). As 
in Hogan (s^ we use scaling to normalize to the high 
frequency results from R.R. Caldwell and B. Allen l5a |: 
10-«(GAi/10-9)i/V^(7a/5)^/^ where p=l. At 
the high frequency limit of our curves this relation is 
used to normalize the parameterization for the number 
of loops created at a given time. 



C. Loop Lengths 

Loop sizes are approximated by the "one-scale" model 
in which the length of created loops scale with the size 
of the Hubble radius. The infinite strings create loops 
with a size that scales with the expansion, creating a 
distribution with a characteristic size a/ H . 

The expression for the average length of newly formed 
loops is {L{t)) = acH{t)~^, where a measures the length 
of the loop as a fraction of the Hubble radius. It is as- 
sumed that the newly formed, stably orbiting loops form 
a range of lengths at a given time, and this range peaks 
near a. For the numerical calculations, and to facilitate 
comparison with previous work, we use a delta function 
for the length. Thus only one length of loop is created at 
a given time and this is the average length defined above. 
Because of the averaging introduced by the expansion, 
this introduces only a small error unless the distribution 
of sizes is large, on the order of several orders of magni- 
tude or more I55j, which would no longer be a one scale 
model. As a check, models of loop formation with sev- 
eral loop sizes are also analyzed, and it is found that the 
larger loops tend to dominate the background over the 
smaller. 

The loops start to decay as soon as they are created 
at time tc, and are described by the equation: 



from Eqn. [T] 



L(tc,t) = acH{tc 



t~tr 



(7) 



(This expression is not valid as L — > but is an adequate 
approximation. As the length approaches zero Eqn.[7]be- 
comes less accurate, but in this limit the loops are small 
and contribute only to the high frequency red noise region 
of the gravitational wave background. This region is very 
insensitive to the exact nature of the radiation so Eqn. [7| 
is accurate within the tolerances of the calculations.) 

Loops that form at a time decay and eventually dis- 
appear at a time td- This occurs when L(tc,td) = and 



= acH{tcy^ -jGfi 



td — tc 



td — tc 



(8a) 
(8b) 



At some time t = td there are no loops created before 
time tc, so they no longer contribute to the background 
of gravitational waves. This is discussed in more detail 
below. 



D. Loop Number Density 

For the number density of strings we parameterize the 
number created within a Hubble volume at a given time 
by Nt/a. The newly formed loops redshift with the ex- 
pansion of the horizon by the inverse cube of the scale 
factor a{t)^^: 



n{t,t') = 



Nt fH{t)Y f a{t) 
a{t') 



(9) 



n(i, t') is the number density at time t as seen by an ob- 
server at time t' , and the redshift has been incorporated 
into the function. In the Appendix there is a descrip- 
tion of the calculation of the scale factor. Analyzing the 
high frequency background equation given by C.J. Hogan 
(32j from R.R. Caldwell and B. Allen [55] the curves are 
normalized by setting Nt — 8.111. 

At an observation time t' the total number density of 
loops is found by summing over previous times t exclud- 
ing loops that have disappeared. 



ns{t') = Vn(t,t') 



(10) 



where tc is the time before which all the loops have de- 
cayed when making an observation at time t' . The earli- 
est time in our sum, tg, is found by solving Eqn. [71 



t' =tc + 



ac^H{tc) 



(11) 



at a given t' . Thus tc{t') tells us that at some time t', 
all of the loops formed before te(t') have decayed and are 
not part of the sum. For large loops and light strings 
tc trends toward earlier in the history of the universe 
for a given t' , allowing more loops to contribute longer 
over cosmic history. For the final calculations, which are 
described in detail in the next section, t' and the number 
density are summed over the age of the universe. 

As the loops decay the energy lost becomes gravita- 
tional radiation, which persists in the universe and red- 
shifts with expansion. This gravitational wave energy 
is detectable as a stochastic background from the large 
number of emitting loops. 
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III. GRAVITATIONAL RADIATION FROM A 
NETWORK OF LOOPS 

A. Frequency of Radiation 

The frequency of radiation is determined by the length 
of the loop, and has normal modes of oscillation denoted 
by n. The frequency of gravitational radiation from a 
single loop is given by, 



2nc 
L{U, t) ' 



(12) 



where tc is the creation time of the loop and t is the 
time of observation. There is also redshifting of the fre- 
quency with expansion which must be accounted for as 
this strongly affects the shaped of the background curves. 
The strings radiate most strongly at the fundamental fre- 
quency so n = 1 for most of the calculations. This also 
facilitates computation times and introduces only small 
errors in the amplitude of the background, which is shown 
in Results. Ultimately the higher modes do not greatly 
alter the results at the resolution of the computations. 



B. Gravitational Wave Energy 

A network of loops at a given time t' radiates gravita- 
tional wave energy per volume at the rate: 



dt 



(13) 



Recall that Ugit') is the number density n{t,t') summed 
over previous times i at a time of observation t'; given 
by Eqn. [TUl 

To find the total gravitational wave energy density at 
the present time the energy density of gravitational waves 
must be found at previous times t' and then summed 
over cosmic history. Since the energy density is given by 
Eqn. [T3] the total density is found by summing over Usit') 
for all times t' while redshifting with the scale factor, a^^. 
Thus the total gravitational energy per volume radiated 
by a network of string loops at the current age of the 
universe is given by: 



ait') 



ns{t')dt', (14a) 



n{t,t') dt'. (14b) 



Here tuniv is the current age of the universe and to is 
the earliest time of loop formation. The internal sum is 
over t and the overall integral is over t' . In practice we 
can make ta small to observe higher frequency, H loops, 
earlier during the radiation era. 

Eqs. [14] do not have explicit frequency dependence, 
which is needed to calculate the background gravita- 
tional wave energy density spectrum Qgwif) determined 



by Eqn. [2l This requires the gravitational energy den- 
sity be found as a function of the frequency: the number 
density as a function of time is converted to a function of 
frequency at the present time. This is discussed in detail 
in the Appendix. 

For most calculations it is assumed each length of loop 
radiates at only one frequency, / = 2c/L. The number 
density, as a function of frequency and time t' , is red- 
shifted and summed from the initial formation time of 
loops to to the present age of the universe: 



Pgwif) = '^Gp'^c 



ait'Y 



Mf- 



ait') 



to 



,t')dt'. 



(15) 

Note that the frequency is also redshifted. 

We have checked the effects of including higher mode 
radiation from each loop. The time rate of change of the 
gravitational energy density of loops at time t' is found 
by summing over all frequencies fj and weighting with 
the power coefficients Pj: 



dpgu,U',t') 

dt' 



Y,PjGp'cn{f;,t'). (16) 



This is then integrated over cosmic time and redshifted 
to give the current energy density in gravitational waves: 



Pgwif) = Gp'^c / dt 



'^4 



Y^P^nif,-^,t')}. (17) 



^ituniv) 



The power coefficients Pj are functions of the mode and 
have the form Pj oc j"'*^'^ [U, [1^. This includes the 
behavior of cusps and kinks, which also contribute an 
/"^/■^ tail to the power spectrum. Depending upon the 
percentage of power in the fundamental mode, the sum 
to infinity is not needed. As an example, if 50% of the 
power is in the fundamental mode only the first six modes 
are allowed. In general it is found that the higher modes 
do not significantly effect the values of the background 
radiation. Details are given in Results. 

From the gravitational radiation density we can find 
the density spectrum readily by taking the derivative 
with respect to the log of the frequency and dividing by 
the critical density: 



^gw if) 



f_dfgJJ)_ 
Pc df 
SnGf dpg^if) 

3iJ2 df • 



(18a) 
(18b) 



A convenient measure is to take ilg^h^ to eliminate un- 
certainties in Ho = Hituniv) — lOO/i km s^^Mpc^^ from 
the final results. 
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IV. RESULTS 
A. General Results 

All of the computations show the expected spectrum 
for loops: old and decayed "H" loops leave a flat spec- 
trum at high frequencies and matter era "P" loops con- 
tribute to a broad peak, with smooth transition region 
connecting them that includes comparable contributions 
from both. The frequency of the peak is strongly depen- 
dent upon the string tension, with lighter strings leading 
to a higher frequency peak. This is a result of the lighter 
loops decaying at a lower rate, so smaller loops survive 
for a Hubble time. This behavior is exhibited in Fig. [TJ 

Fig. 2 shows the peaks of power output per unit vol- 
ume per log frequency as a function of frequency at differ- 
ent times in cosmic history. The most recent times con- 
tribute to the broad peak, while later times contribute 
to the flat red noise region. Although the power output 
per volume is larger earlier, the summation time is fore- 
shortened so the contributions overlap to form the lower 
amplitude flat region in the background. 

The differences in these calculations from previous 
publications can be accounted for by the use of larger 
loops and lighter strings as predicted in [3^ . The change 
in cosmology, i.e. the inclusion of a cosmological con- 
stant, is not a cause of significant variation. Since the 
cosmological constant becomes dominant in recent times 
the increased redshift and change in production rate of 
loops does not strongly affect the background. 



B. Varying Lengths 

The effect of changing the size of the loops is shown 
in Fig. 3. A larger a leads to an overall increase in 
the amplitude of fig^ but a decrease in the amplitude of 
the peak relative to the flat portion of the spectrum. As 
expected there is little if no change in the frequency of the 
peak. In general the red noise portion of the spectrum 
scales with a^l"^ and this is seen in the high frequency 
limits of the curves in Fig. 3 [33]. An increase in a leads 
to longer strings and an increase in the time that a given 
string can radiate. This leads to an overall increase in 
amplitude without an effect on the frequency dependence 
of the spectrum peak. 

In Fig. [4] a smaller a is plotted and shows the reduction 
in the background. The general shape of the background 
remains unchanged, except for some small differences in 
the heavier strings. For heavier strings the amplitude 
of the peak tends to increase relative to the red noise 
portion of the background spectrum, and the frequency 
of the peak shifts less with G/i. These differences become 
more pronounced for heavy and very small strings; and 
at this limit our results are found to match very closely 
those computed in jssj . 



C. Luminosity 

A residual factor of the order of unity arises from un- 
certainty in the typical radiation losses from loops when 
all modes are included averaged over an entire popula- 
tion. Changing the dimensionless luminosity 7 affects the 
curves as shown in Fig. O Larger 7 decreases the ampli- 
tude of the red noise region of the spectrum, but increases 
the lower frequency region. Smaller 7 increases the am- 
plitude of the peak and the flat portion of spectrum, while 
leading to a decrease in the high frequency region. This 
occurs because the loops decay more quickly for larger 
7 and have less time to contribute over the long time 
periods of cosmology. In the matter dominated era the 
higher luminosity allows the loops to contribute more to 
the high frequency region, although it does not increase 
the peak amplitude. Higher luminosity also decreases 
the frequency of the overall peak by contributing more 
power in current loops. The power and change in length 
depend on the luminosity by -E oc 7 and L oc —7. Note 
that the gravitational wave power depends on the square 
of the mass density but the first power of the luminosity, 
thus different behavior is expected varying 7 or fi. An 
analytic description of the scaling of the amplitude and 
frequency of f2g^ with 7 is given in [3^ . 



D. Multiple Lengths 

Some models indicate several characteristic scales for 
the loop populations, with some combination of several 
scales. This implies multiple distributions of lengths, 
which can be described as a combination of different a's. 
Figures [1] [31 and [3] indicate the eventual dominance of 
the larger loops over the smaller in the final spectrum. 
Even though, in this model, the number of loops created 
goes as cx~^, their contribution to the background is still 
much smaller and scales as a^^^. Thus an admixture of 
loop scales tends to be dominated by the largest stable 
radiating loops, which is shown in Fig. [S) 

In Fig. [6] we assume 90% of the loops shatter into 
smaller sizes a — 10~^ while 10% are larger with a = 0.1. 
A "two scale" model is used and two relatively narrow 
distributions are assumed, as in the one scale model. It 
is imagined that the inital loops are large and splinter 
into smaller loops, with the percentage that shatter de- 
pedent upon the model. Overall the amplitude decreases 
by an order of magnitude below a pure a = 0.1 ampli- 
tude, but it is an order of magnitude above the a = 10^^ 
spectra. Thus the large loops, even though they form a 
small fraction, contribute more to the background. Be- 
cause of this weighting a one-scale model is likely to be 
suitable if appropriately normalized. Note that a string 
simulation requires a dynamic range of scale greater than 
a to make a good estimate of a and to choose the right 
normalization. 
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FIG. 1: The gravitational wave energy density per log frequency from cosmic strings is plotted as a function of frequency for 
various values of Gfi, with a = 0.1 and 7 = 50. Note that current millisecond pulsar limits have excluded string tensions 
Gyb > 10~^ and LISA is sensitive to string tensions ~ 10~^^ using the broadband Sagnac technique, shown by the bars just 
below the main LISA sensitivity curve. The dotted line is the Galactic white dwarf binary background (GCWD-WD) from A.J. 
Farmer and E.S. Phinney, and G. Nelemans, et al., [toI. [7ll|. The dashed lines are the optimistic (top) and pessimistic (bottom) 
plots for the extra-Galactic white dwarf binary backgrounds (XGCWD-WD) [713 . Note the GCWD-WD eliminates the low 
frequency Sagnac improvements. With the binary confusion limits included, the limit on detectability of Gfi is estimated to be 



E. High Frequency Modes, Kinks, and Cusps 

Calculations are done which include higher mode be- 
havior of the loops; this sends some of the gravitational 
wave energy into higher frequencies for each loop. Adding 
higher modes tends to smooth out, to a larger degree, the 
curves and slightly lower the peaks on the spectra. When 
applied to the power per volume per log frequency as in 
Fig. [5] we see a less steep high frequency tail at each time 
and a slight decrease in the amplitude of the peak. From 
Fig. [7] a comparison of the fundamental mode versus the 
fundamental mode with higher modes included shows: a) 
the peak decreases, b) the peaked portion widens slightly, 
and c) the red noise amplitude remains relatively un- 
changed. 

Fig. [7] demonstrates a calculation in which ^ 50% of 
the power is in the fundamental mode and the power 
coefficients go as n"''/^ in the higher modes. For the 
fundamental mode Pi = 25.074, and the next five modes 
are P2 = 9.95, P3 = 5.795, P4 = 3.95, P5 = 2.93, and 
Pe = 2.30. Their sum is -P™ = 7 = 50. The amplitude 
of the peak shifts from 1.202 x 10"^ to 1.122 x 10"^ a 
6.7% decrease. The frequency of the peak shifts more 



substantially, from 3.162 x lO^'^ Hz to 7.080 x 10^"^ Hz, 
an increase of 124%. This increase in frequency of the 
peak does not affect the limits of detection by LISA be- 
cause the amplitude shrinks very little. Note that the 
high frequency red noise region is essentially unchanged, 
and for Gfi = 10~^^ this is the region in which LISA is 
sensitive. 

Plots of spectra for other G/i with high frequency 
modes show changes similar to those in Fig. [T] The peak 
amplitude drops slightly and the frequency of the peak 
is increased. Overall, the limits from Table HTl remain un- 
changed by the introduction of 50% of the power in the 
higher frequency modes. 

Cusps and kinks in the loops are responsible for high 
frequency "bursts" of gravitational wave energy [47l , |48| . 
l49l. Isoj. They are two different manifestations of string 
behavior: cusps are catastrophes that occur approxi- 
mately once every oscillation, emitting a directed burst 
of gravitational radiation; while kinks are small wrin- 
kles which propagate on the loops, emitting a higher fre- 
quency of directed gravitational energy. The light strings 
make it less likely to detect these random bursts; in par- 
ticular they are more difficult to pick out from their own 



8 



|y=50, ti'=0.1, N,=8.111,G^=10-'^| 



(y=50, ff=10-\N,=8.111} 



10 



E 



o I „ 









' \__t=10"' yr 










-/X/ t=13.7Gyr 





-10 



-6 -4 -2 
log (freq in Hz) 



5 
G 




gcwd-;!vd/ /: xgcwd-wd 

-6 -4 -2 

log (freq Hz) 



FIG. 2: The power per unit volume of gravitational wave 
energy per log frequency from cosmic strings is measured at 
different times in cosmic history as a function of frequency 
for G/i — 10~^^, a — 0.1, and 7 — 50. Note that early in the 
universe the values were much larger, but the time radiating 
was short. The low frequency region is from current loops 
and the high frequency region is from loops that have decayed. 
This plot shows the origin of the peak in the spectrum at each 
epoch; the peaks from earlier epochs have smeared together 
to give the flat high frequency tail observed today. 
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The gravitational wave energy density per log fre- 
from cosmic strings is shown with varying a for 



FIG. 3: 

quency 

G/i — 10~^^ and 7 — 50. Here a is given values of 0.1, 10"'^, 
and 10~® from top to bottom. Note the overall decrease in 
magnitude, but a slight increase in the relative height of the 
peak as a decreases; while the frequency remains unchanged. 
The density spectrum scales as ci'^'^ , and larger loops domi- 
nate the spectrum over small loops. 



confusion background [32|- Since cusps, kinks, and higher 
modes add an tail to the power spectrum of each 

loop, it is reasonable to expect that they will have a simi- 
lar effect on ■ The general behavior of high frequency 
contributions is represented in Fig. [T] 



FIG. 4: The gravitational wave energy density per log fre- 
quency from cosmic strings for a — 10~® and 7 = 50. Smaller 
a reduces the background for equivalent string tensions, with 
a LISA sensitivity limit of G/i > 10~^* if the Sagnac sensi- 
tivity limits are used. The bars below the main LISA sen- 
sitivity curve are the Sagnac improvements to the sensitiv- 
ity, while the dotted line is the Galactic white dwarf binary 
background (GCWD-WD) from A.J. Farmer and E.S. Phin- 
ney, and G. Nelemans, et al., [toI. [71||. The dashed lines are 
the optimistic (top) and pessimistic (bottom) plots for the 
extra- Galactic white dwarf binary backgrounds (XGCWD- 
WD) [7l[. Note the GCWD-WD eliminates the low frequency 
Sagnac improvements, and increases the minimum detectable 
G/i to > 10^^^. G/i ranges from the top curve with a value 
10"'^ to the bottom curve with lO"^'^. 



V. MILLISECOND PULSAR LIMITS 

A. Limits on Background 

Millisecond pulsars have been used as sources to indi- 
rectly measure a gravitational radiation background for 
a number of years [H, [H, [H, [13, [H, HI] . Limits on the 
strain, and correspondingly the energy density spectrum 

at three frequencies: 
|. The limits on the 
strain vary with the frequency dependence assumed for 
the source of gravitational radiation at the corresponding 
frequencies: 



figuj, have been estimated recently 
l/20?/r-\ l/8 2/r-\ and 1 yr"^ %. 



hcif) - A 



f 



yr- 



13 



(19) 



where (3 is the power of the frequency dependence of the 
source and is denoted as a in F.A. Jenet, et al., 2J]. A 
table of values is provided in that source for various cos- 
mological sources of stochastic gravitational radiation, 
with /? — —7/6 used for cosmic strings. From the 
cosmic string limits are: 



gw 



1 yr-')h' 



f]g„(l/20 yr-^)h 



< 9.6 X 10-^ (20a) 
^)h'^ < 1.9 X 10"^ (20b) 
2 < 2.6 X 10~^ (20c) 
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FIG. 5: The energy density in gravitational waves per log 
frequency from cosmic strings with varying luminosity 7 for 
G/i = lO"'^^ and a = 0.1. Here the 7's are given as follows 
at high frequencies: 50 for top curve, 75 for the middle curve, 
and 100 for the bottom curve. The peak frequency, peak 
amplitude, and amplitude of flat spectrum all vary with the 
dimensionless gravitational wave luminosity 7. 



|7=50, Q'=(10% 0.1 & 90% 10-'), N,=8.111) 

LISA 




-8 -6 

log (freq Hz) 

FIG. 6: The gravitational wave energy density per log fre- 
quency from cosmic strings for 7=50, G/i = 10"^^, and a 
combination of two lengths (q's). 90% of the loops shatter 
into loops of size a = 10~^ and 10% remain at the larger 
size a = 0.1. The overall decrease is an order of magnitude 
compared to just large alpha alone. From bottom to top the 
values of G/i run from 10"" to 10"*. 



Results of our calculations show that for heavier strings 
at the measured frequencies /3 is approximately -1, but 
for lighter strings it varies from this value. The variation 
of /3 is summarized in Table U which has been calculated 
at j — 1/20 yr^^ and f — 1 yr^^ ■ This change in (3 will 
result in limits differing slightly only at low frequencies 
from those in Eqs. HOI [24]. 

In Figures [T] and 2] the limits for millisecond pulsars 
[23 | are shown, with allowed values of f2g^ indicated by 
the direction of the arrow. For large strings a = 0.1, the 
maximum allowed Gji is about 10"^, Fig. [TJ Decreasing 
a decreases the amplitude of the background allowing for 




-6 -4 

log (freq Hz) 

FIG. 7: The gravitational wave energy density per log fre- 
quency for 7=50, G/x = 10"^'^, a = 0.1; plotted in one case 
(dashed) with only the fundamental mode, while for another 
case (solid) with the first six modes. The dashed curve has 
all power radiated in the fundamental mode, while the solid 
curve combines the fundamental with the next five modes: 
Pi = 25.074, P2 = 9.95, P3 = 5.795, P4 = 3.95, P5 = 
2.93, Pq = 2.30. Note the slight decrease in the peak and its 
shift to a higher frequency. The LISA sensitivity is added as 
a reference. 



heavier strings within the pulsar measurement limits. A 
more complete list is displayed in Table |TT1 



B. Dependence of Characteristic Strain on 
Frequency 

From the curves of Fig.[T]one can find the frequency de- 
pendence of the characteristic strain as given by hc{f) = 
A{f I yr~^Y . Again, (3 is the power law dependence of 
he on frequency, and is generally written as a with a nu- 
merical value of -7/6 [2J, [63, [IJl . In this paper (3 is used 

" is 



to prevent confusion with the length of loops, 
known, then using Eqn. [23]wc find: 



If n 



gw 



^giu (/) 



^2(1- 



3iJ2 (3.17 X 10- 

/3 



(21) 



With this equality (3 is found by computing 
d\n{ngy,)/d\nf at / = 1 yr-^ = 3.17 x IQ-^Hz 
and setting this equal to 2(1 + /3). The results are listed 
in Table U 

From Fig. [1] the change in slope with Gfi is evident at 
the measured frequency /o = 1 yr~^. For large values of 
G^, fo is in the flat red noise section of the spectrum so 
the slope of ilgw goes to zero, and f3 approaches -1. As 
the string tension is reduced, (3 is measured first at the 
high frequency tail so the slope is negative, and thus /3 is 
more negative. Once over the peak on the high frequency 
side, the slope becomes positive, thus /3 must be greater 
than -1. The characteristic strain power law dependence 
(3 then approaches a value of -0.67. 
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TABLE I: P for given Gfi at f = 1/20 yr'^ and / = 1 yr'^ , 
a = 0.1, 7 = 50, where he oc f'^. 





/ = 1/20 yr-i 


/ = lyr-i 


10"^ 


-1.00306 


-1.00031 


10"® 


-1.05314 


-1.00894 


10-^ 


-1.05985 


-1.03853 


10-10 


-0.931749 


-1.06593 


10-" 


-0.735778 


-0.987882 


10-12 


-0.682874 


-0.778679 


10-13 


-0.676666 


-0.688228 


10-14 


-0.676034 


-0.676054 


10-1^ 


-0.675971 


-0.674794 


10-16 


-0.675965 


-0.674668 



This dependence on G/i changes the value of the con- 
stant A in Eqn. [inland thus the hmits on [23|- Higher 
mode dependence only modestly changes the values of (3 
since it approaches the limits -1 and -2/3 for high and 
low frequencies respectively. The frequency at which it 
takes on these values depends upon G/i as seen in Fig. [1] 

VI. FUTURE LISA SENSITIVITY 
A. Physical Principles 



TABLE 11: Millisecond Pulsar Limits and LISA sensitivity for 
7 = 50. WMAP and SDSS 35] have excluded Gfi > 10~\ 



a 


Millisecond Pulsar 
Limits, Gfi < 


Minimum LibA 
Sensitivity, Gfi > 


0.1 


w-'> 


10-16 


10"2 


io-« 


10-15 


10-3 




10-14 


10-* 




10-13 


10-^ 




10-12 


10-6 




10-11 



^gw{f), with the strain that LISA will detect. This limit 
on the detectable strain is shown as a curve on our plots 
and is taken from (67j . 

In this paper the LISA sensitivity curve is cal- 
culated with the following parameters: SNR=1.0, 
arm length=5xl0^ m, optics diameter=0.3 m, laser 
wavelength=1064 nm, laser power=1.0 Vl^, optical train 
efRciency=0.3, acceleration noise=3 x 10-i^m/ (s'^^/Wz), 
position noise=2 x IQ-^^ m/^/Hz. 

Extracting the background signal by combining signals 
over a broad band and an extended time can increase the 
sensitivity of the detector. One method known as the 
"Sagnac technique" is shown on plots of the background 
[68ll69j. for example see Fig. [1] 



The pending launch of LISA will give new oppor- 
tunities to observe gravitational waves in general, and 
a stochastic background in particular [2^, l26l [39l . [6^ . 
There are a number of theoretical sources for this back- 
ground, including cosmic strings [13, [H, HE HI]- 
LISA sensitivity has been modeled and methods for im- 
proving sensitivity by calibrating noise and integrating 
over a broad band have been proposed [b^, HI, [6§] . The 
cosmic string background density spectrum calculated 
here is compared to the potential for detection by LISA. 

Detectors such as LISA measure the strain on space- 
time, so a relation between the background spectrum and 
strain is needed. The rms strain is denoted by hrmsif) or 
h{f), and for an isotropic, unpolarized, stationary back- 
ground [2l[3i,[6ll,[6|: 

47r2 

^sUf)^^j'hl„,M)- (22) 

Another measure is the characteristic strain spectrum 
hc{f)^ = 2//i^„g(/), when substituted gives: 

- J^fhlif). (23) 

In the literature, Sfi{f) — h„ns(.f )'^ is often used and 
referred to as the spectral density [61[. These relations 
are important as they relate the energy density spectrum 



B. Limits for LISA Sensitivity 

For both Figures [1] and 2] there are limits on the de- 
tectability of the background due to strings of various 
mass densities by LISA from [63|. From Fig. [1] with 
a = 0.1 the minimum value of the string tension measur- 
able is ~ IQ--'^^. Fig. 4, which plots a smaller a = 10-^, 
the minimum is G/i ~ IQ-^^. These values and others are 
listed in Table |IT1 Sagnac observables increase our avail- 
able spectrum and change the detectable string tensions 
as indicated [ssj . 

Shown in Table |TT] are the limits on string tensions for 
string sizes a of 0.1 to 10"®, along with the correspond- 
ing minimum tension that LISA can detect. Galactic and 
extra-Galactic binaries confusion limit is taken into ac- 
count for this table. For a < 10~^ the maximum string 
tension will be greater than 10-^, although these val- 
ues are excluded by WMAP and SDSS [1^. It is clear 
from the table and figures that decreasing a increases the 
maximum string tension allowed as well as increasing the 
minimum string tension detectable by LISA. 



C. Confusion Noise from Binary Systems 

Galactic white dwarf binaries add a background "con- 
fusion" that shrouds portions of LISA's spectrum. This 
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background has been studied and predictions given for 
its spectrum [70| . There are also extra-Galactic binary 
sources which similarly create confusion noise [7l| . These 
confusion limits are shown on plots of the gravitational 
radiation spectra, and they affect the minimum values of 
string tension that are detectable. 



VII. CONCLUSIONS 

Broadly speaking, the results here confirm wit h g reater 
precision and reliability the estimates given in [32] ■ Our 
survey of spectra provides a set of observational targets 
for current and future projects able to detect stochastic 
gravitational wave backgrounds. For light strings, the 
string tension significantly affects the spectrum of gravi- 
tational radiation from cosmic string loops. In addition 
to a reduction in the overall radiation density, the most 
significant change is the shift of the peak of the spectra to 
higher frequencies for light strings. Near the current pul- 
sar limits, the peak happens to lie close to the frequency 
range probed by the pulsars, so spectra as computed here 
are needed as an input prior for establishing the value of 
the limits themselves. For much lighter strings, the peak 
is at too high a frequency for strings to be detectable by 
pulsar timing as the spectrum falls below other confusion 
backgrounds. Lighter strings will be detectable by LISA. 
For the lightest detectable strings the peak happens to lie 
in the LISA band and the detailed spectrum must again 
be included in mounting observational tests. The spec- 
trum from strings is quite distinctive and not similar to 
any other known source. 

A high probability of forming loops with stable, non- 
self-intersecting orbits, as suggested by recent simula- 
tions, leads to larger string loops at formation which 
in turn give an increased output of gravitational wave 
energy and improved possibilities for detection for very 
light strings. Recent simulations have shown stable radi- 
ating loops form at a significant fraction of the horizon, 
of order a = 0.1; for loops of this size our calculations 
show that the maximum string tension allowed by current 
millisecond pulsar measurements is G/i < 10~^, and the 
minimum value detectable by LISA above estimated con- 
fusion noise is Gfi « 10~^^. In field theories, the string 
tension is typically related to the scale by G/z oc A^ / m^; , 
so the maximum detectable scale in Planck masses cur- 
rently allowed by millisecond pulsars is A < 10~^'^, or 
around the Grand Unification scale; with LISA, the limit 
wiU be about A^ 10"* or 10" GeV. These results 
suggest that gravitational wave backgrounds are already 
a uniquely deep probe into the phenomenology of Grand 
Unification and string theory cosmology, and will become 
more powerful in the future. The most important step in 
establishing these arguments for an unambiguous limit 
on fundamental theories will be a reliable quantitative 
calculation of a. 



APPENDIX: NUMERICAL CALCULATIONS 
1. General 

For computations the relevant differential and integral 
equations are solved numerically. Here is a list of the 
equations that are used, most of which have been dis- 
cussed in the sections above. In addition there are several 
more that are expanded upon in this section. Below note 
that n{t,t') is the number density with time t at an ob- 
servation time t' , while n{f',t') is the number density as 
a function of frequency and time t' . The time t is always 
earlier than the time t' , which is the time of observation 
of the quantity in question. 



SttG 



(A.l) 



P — Prad + Pi/ + Pmatter + pA 
H~ Poo ^~ PL H" Pgw ; 

(A.2) 



fM = 



2c 



Lit,t'y 



(A.3) 
(A.4) 



(A.6) 



a \ c J \a{t') 



Pgwif) = iGfi C I 



X n{f 



ait') 



,t')dt' 



Recent observations allow us to neglect loops pL , infinite 
strings poo , and gravitational waves pgw in the Friedmann 
equation, so these terms are removed (72i] . 



2. Scale Factor 

Calculations of the scale factor are done using the stan- 
dard cosmological model with a homogeneous isotropic 
spacetime defined by the Friedmann-Robertson- Walker 
metric, solved using Einstein's equations [Ijlll, [ill- As 
an addition the effects of a cosmological constant on the 
cosmic string gravitational wave background are com- 
puted. The energy densities used to describe the dynam- 
ics of the scale factor are radiation, neutrinos, matter, 
and a cosmological constant. The equation of motion for 
the scale factor is determined by the Friedmann equation. 



SttG \ ^ / \ r,\ 

—^{Prel + Pm) + , (A.7) 



where the p's are functions of the scale factor. We then 
integrate the above equation to find a{t). 

The relativistic and matter densities follow the rela- 
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tions: 



Prel 
Prn — 



{Prel)oat 
{Pni)oal 



(A.8a) 
(A.8b) 



The gravitational radiation density also goes as a(t)~^. 

In order to compute the scale factor the Friedmann 
equation is transformed into an integral equation and 
solved numerically. The main source of error is the rela- 
tivistic treatment of neutrinos for all time, but this will 
incur a very small error in more recent times. Given that 
neutrinos have mass, it is possible several species will be 
non-relativistic at the present time. This is a minimal 
consideration as the matter and dark energy contribu- 
tions are larger by at least five orders of magnitude when 
the neutrinos are non-relativistic. The following param- 
eters are taken from the Particle Data Group, 2006 [73 |: 





= 0.76, 


(A.9a) 


^'matter 


= 0.24, 


(A.9b) 


^rad 


= 4.6 X 10-^ 


(A.9c) 




= 3.15 X 10-^ 


(A.Qd) 


h 


= 0.73, 


(A.9e) 


Pc 


SttG" 


(A.9f) 



where Vlj — Pj/Pc, and the constant h is given by Ho 
lOO/i km s-^Mpc'^. 



3. Gravitational Radiation from Loops 

At a given time t', there is an earlier time t before 
which all of the loops have decayed, which we denote t^- 
To find this value at any time t' , Eqn. [7] is set to zero 
giving L(t,t') — 0, and t — t,, is calculated numerically. 
Note that is a function of t'. For very early times te{t') 
can be lower than the initial time in our sum over cosmic 
history Iq. In this case to was used instead of the value 

for te{t'). 

For many of the calculations it is assumed each loop 
radiates at only one frequency, / — 2c/ L. Then a nu- 
merical list of values is computed for t — to t — t' , 
{L{t, t'), n{t, t')}, and an interpolation function is created 
from the list. This is the number density as a function of 
length and time t' , n{L, t'). Next a list of frequency and 
number density is created, {2c/L,n{L,t')}, with lengths 
running from L(te, t') to L(t' , t'). An interpolation func- 
tion is made from this list, and this results in the number 
density as a function of time t' , n{f' , t'). But the frequen- 
cies must be redshifted to the current time tuniv so the 
function is n{fa(t')/a{tuniv),t'). 

Then we must integrate over all t' from the earliest 
time of loop formation to the present age of the uni- 



verse tuniv This integral. 



a{t') 



is transformed into a sum, 



(A.IO) 



Pgw{f) = iGp 



t'=t„ 



aitr 



nif- 



ait') 



-,t')Ai'. 



(A.ll) 

Suitable At"s are chosen to ensure both computability 
and to minimize errors. In general we solve all equations 
in logarithmic bins so At' will vary in accordance. Using 
small values of At to facilitate computation speed will 
result in relatively small error; choosing At = 0.1 in log 
scale gives a 10% smaller value for ftg^ compared to At = 
0.001 or At = 0.01 in log scale. Values smaller than 
At = 0.01 give little or no change to the background. 
After the sum the result is the gravitational wave density 
as a function of frequency at the present time. 

To include higher modes the power coefficients Pj must 
be included in the calculation. The subscript j is used 
instead of n here to prevent confusion with the number 
density, n(/',t'). At each time t' the number density, 
and thus the time rate of change of the energy density, of 
loops is calculated exactly as above, but more frequencies 
are included. Each number density at frequency fj is 
summed and weighted with Pj , and this number density 
is then integrated over cosmic history, 

p,^if) = Gp'c dt'^i^ 

J to ^{J'univ ) 

x<E^."(/.7#^'0^ (A.12) 



The overall integral is then transformed into a sum as 
before. The sum of the power coefficients is Pj = 7i 
where the Pj (x 

The final result is to calculate flgwif) at the present 
time. This is done by finding the current energy density 
of gravitational waves as a function of frequency, and 
then taking the derivative with respect to f: 



^gw (/) 



/ dp 



gw 



Pc df 



(A.13) 



This is then plotted as log{ilgyj{f) h?) vs. logif /Hz), 
as seen in Figures 1,3—7. This quantity is, again, the 
density of stochastic gravitational wave energy per log of 
frequency in units of the critical density and is related to 
the strain on space by the gravity waves. 
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